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ABSTRACT 

We discuss V and R band photometry for 67% of the SuUivan et al. (2000) SA57 ultraviolet-selected 
galaxy sample. In a sample of 176 UV-selected galaxies, Sullivan et al. (2000) find that 24% have (UV-i3) 
colors too blue for consistency with starburst spectral synthesis models. We propose that these extreme 
blue, UV excess galaxies are Wolf-Rayet (WR) galaxies, starburst galaxies with strong UV emission 
from WR stars. We measure a median {V — R) = 0.38 ± 0.06 for the UV-selected sample, bluer than a 
sample optically selected at R but consistent with starburst and WR galaxy colors. We demonstrate that 
redshifted WR emission lines can double or triple the flux through the UV bandpass at high redshifts. 
Thus the (UV— -B) color of a WR galaxy can be up to 1.3 mag bluer at high redshift, and the expected 
selection function is skewed to larger redshifts. The redshift distribution of the extreme blue, UV excess 
galaxies matches the selection function we predict from the properties of WR galaxies. 

Subject headings: galaxies: starburst, galaxies: photometry, stars: Wolf-Rayet, ultraviolet: galaxies 



1. INTRODUCTION 

Sullivan et al. (2000, hereafter S2000) published an ul- 
traviolet (UV) flux-limited catalog of galaxies, imaged 
with the FOCA balloon-borne UV camera of Milliard et al. 
(1992). The catalog is an extension of work by Treyer et al. 
(1998), and contains 273 galaxies with redshifts z < 0.4. 
Vacuum UV imaging enables the direct detection of radi- 
ation from hot, massive stars and, by association, actively 
star forming galaxies. The luminosity function (LF) of the 
UV-selected galaxies is somewhat surprising. S2000 find 
a faint end slope in the Schechter (1976) parametrization 
of the LF, auv = —1-55 ± 0.11, steeper than typical opti- 
cal LF faint end slopes (Folkes et al. 1999; Marzke et al. 
1998; Geller et al. 1997, etc.). As a resuh, S2000 find no 
strong decrease in UV luminosity density between 2 ~ 0.4 
and z = and thus little decline in the density of star 
formation over this range of cosmic epochs. 

In contrast, Madau et al. (1996) conclude that there is a 
factor of ~ 10 decline in the comoving density of star for- 
mation between redshifts z = 1 and z = 0. This conclusion 
is based on a combination of rest frame near-UV luminosi- 
ties from the /-selected Canada France Redshift Survey 
(Lilly et al. 1995) and local Ha measures from Gallego et 
al.'s (1995) objective prism survey. 

Curiously, a substantial fraction of the S2000 UV- 
selected galaxies have (UV—B) colors bluer than the bluest 
color that can be produced by starburst spectral synthesis 
models (Fioc & Rocca-Volmerange 1997; Leitherer et al. 
1999). S2000 could therefore not explain the color of these 
galaxies. 

Here we propose that the extremely blue galaxies are 
Wolf-Rayet (WR) galaxies. WR galaxies show direct sig- 
natures of WR stars in their spectra, usually a broad He 
II A4686 feature originating in the stellar winds of WR 
stars. This subtle feature is not seen in the fiber-spectra 
of S2000 (M. Sullivan 2000, private communication), but 
requires good S/N and is spatially dependent on the in- 
clusion of starburst regions (Schaerer et al. 1999). The 
starburst regions dominate the broadband UV fiux and. 



though short-lived (< 10^ yrs), WR stars can dominate 
the UV light of young starbursts. 

In this paper we describe optical photometry for 67% 
of the S2000 Selected Area 57 (SA57) UV-selected galaxy 
sample. Our goal is to compare the optical and UV colors, 
and to suggest that WR emission lines are the source of 
the extremely blue colors found by S2000. In section 2 we 
discuss our data, in section 3 we discuss the (V — R) colors, 
and in section 4 we discuss the impact of Wolf-Rayet UV 
emission lines on the interpretation of the UV sample. We 
conclude in section 5. 

2. DATA 

We obtained broadband V and R images in 1998 Decem- 
ber and 1999 February with the MOSAIC camera on the 
Kitt Peak National Observatory 0.9m telescope. MOSAIC 
is an 8 CCD array with pixel scale 0.424 arcsec pixel"^ and 
field of view 59' x 59' on the 0.9m. We obtained three off- 
set 200 sec R and 333 sec V images at each position along 
a strip in RA. Our data serendipitously overlap with 148 
of the 222 S2000 UV-selected galaxies in the SA57 region. 

We processed the images in the standard way with the 
IRAF MSCRED package. Total magnitudes are obtained 
using SExtractor (Bertin &Arnouts 1996). Seven percent 
of the magnitudes are obtained with a modified version 
of IRAF-based GALPHOT (Freudling 1993). GALPHOT 
allows photometry of galaxies poorly fit by SExtractor, 
including low surface brightness galaxies or galaxies with 
close companions. A comparison between the two software 
packages yields a 1-a scatter in the total galaxy magni- 
tudes of ±0.05 mag. We use Landolt fields (Landolt 1992) 
and the M67 cluster (Montgomery, Marschall, & Janes 
1993) for flux calibration. Depending on the nightly see- 
ing conditions, the error in the standard star calibration is 
±0.02 — 0.03 mag. A direct measure of our intrinsic error 
is the RMS scatter of measurements of identical galaxies 
taken on different nights: ±0.04 mag. 

The photometry includes galactic extinction corrections 
from Schlegel, Finkbeiner, & Davis (1998), who combine 
CORE /DIRBE and IRAS/ISSA data to map dust emis- 
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sion in the Milky Way. Typical extinction values are 
An = 0.025 and Ay = 0.03 mag. The photometry also 
includes "k corrections" for the changing rest frame band- 
pass (Frei & Gunn 1994). Because we know the iV—R) 
color and the redshift of each galaxy, we can determine 
accurate k corrections. Typical values for the V and R 
bandpasses are —0.1 < Ur < 0.1 and —0.15 < ky < 0.15 
mag. 

2.1. Existing UV and B measurements 

The S2000 UV measurements were made from a balloon- 
borne, 40 cm Cassegrain telescope stabilized to 2 arcsec 
RMS. The filter response approximates a Gaussian cen- 
tered at 2015 A, with a FWHM of 188 A. The limiting 
magnitude in the UV catalog is mf/y = 18.5. S2000 state 
that the UV zero-point is accurate to < 0.2 mag; the un- 
certainty in the relative photometry may reach ~ 0.5 mag 
at the limiting magnitude. 

With the UV images in hand, S2000 searched for opti- 
cal counterparts in digitized POSS plates. S2000 obtained 
accurate positions and B photometry, with uncertainty 
±0.20 mag, from the POSS plates. We follow S2000 and 
only use galaxies with a single optical counter-part when 
discussing UV colors. S2000 used the resulting (UV— B) 
colors to calculate k corrections. Finally, S2000 apply a 
reddening correction based on Hq and H/'? line strengths 
and the empirical Calzetti (1997) extinction law for star- 
burst galaxies. They applied a mean reddening correction 
Ay = 0.97 to galaxies lacking Ha and/or H/3 emission. 

The (UV-S)o colors we quote from S2000 include k 
corrections and reddening corrections. Our [V — R) col- 
ors include k corrections but no reddening corrections, be- 
cause S2000 do not provide enough information to repro- 
duce their reddening corrections. Even if this reddening 
data were provided, 72% of the galaxies with {V — R) col- 
ors lack both Ha and H/3 emission. For these galaxies, the 
mean reddening correction would add a uniform ~ —0.2 
mag offset to the [V — R) colors. 

3. OPTICAL COLORS 

Figure 1 shows that the [V — R) colors of the UV- 
selected galaxies are significantly bluer than those from 
an i?-selected galaxy sample, the Century Survey (Geller 
et al. 1997). The Century Survey is a complete red- 
shift survey of 102 square degrees to a limiting magnitude 
R = 16.13, containing 1762 galaxies. The Century Sur- 
vey galaxies have median (y — R)cs = 0.58 ± 0.04; the 
UV-selected galaxies have median (V — i?) = 0.38 ± 0.06. 

The distribution in Figure 1 appears double-peaked, 
but application of the non-parametric estimator of Pisani 
(1993) shows the second peak is not significant. The Pisani 
(1993) kernel probability estimator allows estimation of 
the probability density distribution underlying a data sam- 
ple. In applying this we assume an average I-ct {V — R) 
error ±0.066 mag. 

The S2000 (UV— i?)o colors provide a measure of the 
star formation activity, but there is no correlation between 
the {V-R) and (UV-B)o colors. 

Unlike the extreme blue (UV— i3)o colors found by 
S2000, however, the optical {V — R) colors are consistent 
with starburst galaxy colors. For example, the template 
star burst-galaxy spectral energy distributions SBl and 
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Fig. 1. — Histogram of {V — R) colors for UV-selected galaxies 
(solid line) and i?-sclcctcd Century Survey (dashed line). 

SB2 of Kinney et al. (1996) have synthetic {V-R)= 0.26. 
All but two of the S2000 UV-selected galaxies are redder 
than {V — R) = 0.23, and we note that these blue out- 
liers appear atypically faint and point-like in our images. 
Thus ^ 95% of the UV galaxy {V — R) colors arc consistent 
with starburst galaxy colors. If the extreme blue galaxies 
were AGN, their average color should be iV—R) = 0.60 
(Kotilainen et al. 1993), which is clearly not the case. 

4. DISCUSSION 

An unresolved issue from the S2000 paper is the exis- 
tence of galaxies with (UV— i3)o colors < —2, colors too 
blue to be explained by standard spectral synthesis mod- 
els. In our subsample of the S2000 UV-selected galaxies, 
48 of 148 galaxies have (UV-B)o < -2, yet the opti- 
cal {V — R) colors show no significant departure from the 
rest of the sample. A KS test gives an 80% probability 
of the extreme (UV— B)o < —2 and the other UV-selected 
galaxies being drawn from the same distribution of {V — R) 
colors. 

Next wc examine the redshift distribution of the ex- 
tremely blue galaxies. The lower panel of Figure 2 
shows the redshift distribution of the extreme blue 
(UV— i3)o < —2 galaxies, and the upper panel the re- 
maining (UV— i?)o > —2 galaxies. Figure 2 shows that 
the extreme blue galaxies are preferentially found at high 
redshift. A KS-test yields a 10~^^ probability of the two 
populations being drawn from the same distribution. 

The peaks in the redshift histograms (Figure 2) arc well 
matched to the z ^ 0.3 limit of the Medium Deep Sur- 
vey (WillmcT ct al. 1996), which includes the SA57 region. 
The Medium-Deep Survey is ~ 70% complete to Bj ~ 20.5 
mag, and gives a sense of the existing large scale structure 
in the direction of UV survey. Thus the extremely blue 
galaxies appear to trace the same large-scale structure as 
"normal" galaxies. 

4.1. Wolf-Rayet galaxies. 

The very blue (UV— i?)o < —2 galaxies are an extreme 
population of galaxies found mostly at larger redshift. An 
explanation probably requires emission lines redshifting 
into the UV bandpass. Because S2000 and our (V — R) 
color distribution rule out AGN, it seems plausible the 
emission lines come from a stellar emission line source, 
i.e. Wolf-Rayet (WR) stars, planetary nebulae, or sym- 
biotic stars. Because there are known galaxies with WR 
features, we suggest that the extreme blue (UV— B)o < —2 
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Fig. 2. — Redshift histograms for the (UV-S)o > -2 galax- 
ies (upper panel) and the extreme blue (UV— S)o < —2 galaxies 
(lower panel). Note the excess of extreme blue galaxies at large red- 
shift. The selection function curves are calculated with the S2000 
LF parameters (dashed-llne) and the SSRS2 Sm-Irr LF parameters 
(dash-dot, solid lines) as described in section 4.1. 

galaxies arc WR galaxies. The strongest UV features come 
from ionized C, O, He, Si, and Fe in the stellar winds of 
WR stars (e.g. HilHer & Miller 1999). 

Our optical spectra of 10 extreme blue galaxies all show 
narrow emission lines indicative of active star formation. 
WR galaxies are characterized by the direct signatures of 
WR stars in their spectra, but our S/N is too poor to 
detect a possible A4686 WR feature convincingly. 

To calculate a typical (V — R) color for WR galaxies, 
we searched the NASA/IPAC Extragalactic Database for 
photometery of known WR galaxies. The average optical 
color of 13 WR galaxies is (V-R) = 0.27 ±0.34 mag. This 
average color is in agreement with the blue color we find 
for the (UV-S)o < -2 sample, (V-R) = 0.35±0.11 mag. 

Figure 3 illustrates how WR emission lines are redshifted 
into the S2000 UV bandpass. We plot a sample WR ul- 
traviolet spectrum, WC star HD 165763 (Hillier & Miller 
1999), at redshift z = and at redshift z = 0.285 with the 
S2000 UV and B bandpasses overlaid. Integration of the 
UV bandpass over the redshifted spectrum shows that the 
WR emission lines can boost the UV flux by a factor of 
< 3 compared to the flux at z = 0. We also integrate the 
B bandpass over the Hillier & Miller (1999) WC stellar 
spectrum, and find the (UV— J5) color is shifted up to 1.3 
mag bluewards over the redshift range z=0 to z=0.5 (see 
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Fig. 3.— WR spectrum of star HD 165763 at z = and at 
z = 0.285, with rest-frame UV and B bandpasses. 



Figure 4), blue enough to explain 2/3 of the extreme blue 
(UV— 5)0 galaxies in the SA57 UV-selected sample. 

The WC stellar spectrum represents a maximal effect; 
known WR galaxies have typical WC/WN ratios of 0.2 - 
0.4 (Schaerer et al. 1999). We perform the (UV-B) cal- 
culation on the WN stellar spectrum, WR152 (Hamann 
& Koesterke 1998), and find the (UV-B) color shifts up 
to 1.1 mag bluewards over the same redshift range (see 
Figure 4). 

We acknowledge that strong UV emission lines are not 
common features of known low redshift WR galaxies (Lei- 
therer 2000, private communication). We calculate the 
shifts in (UV-B) for a dozen of the McQuade et al. (1995) 
"UV to optical spectral distributions of northern star- 
forming galaxies," and find an average A{\JV — B) ^ —0.3 
mag (mostly due to the continua, not the emission lines). 
Yet good examples of the very blue objects may not be 
found locally. The UV selection should pick out extreme 
objects; known WR galaxies are optically selected. We 
note the very blue UV galaxies have median z=0.2. By 
comparison, the Schaerer et al. (1999b) catalog of known 
WR galaxies drops off sharply at z=0.02. 

We can also take advantage of HST FOS observations of 
the Small Magellanic Cloud H II region N88A (Kurt et al. 
1999) as a UV spectral templates for WR galaxies. When 
the 82000 bandpass is centered on the strong ClII and Silll 
lines in this object, the UV flux doubles compared to the 
surrounding continuum level. 

S2000 compute their luminosity function (LF) with the 
full UV sample; here we construct a model LF for the WR 
galaxies alone. The impact of the emission lines on the 
observed UV flux changes with redshift. We explore the 
selection function for finding WR galaxies in a magnitude 
limited redshift survey by computing the change in abso- 
lute magnitude resulting from the UV emission lines. 

WR galaxies are blue and show active star formation. 
We thus use the Second Southern Sky Redshift Survey 
(Marzke et al. 1998, SSRS2) LF for the bluest galaxies, the 
Magellanic spiral and irregulars, as the basis of our model. 
The SSRS2 is a B-selected survey complete to mi, < 15.5 
that contains 5404 galaxies, all with morphological clas- 
sifications. The SSRS2 Sm-Irr LF has a steep faint-end 
slope, a = -1.81, similar to the S2000 LF, a = -1.55. 

We compute the selection function for the WR galaxies 
according to equation (1), 

N{< Miirr,) oc / Dc{zfdDc{z) / cj>{M,z)dM 

(1) 

where Dc{z) is the co-moving distance. We assume a fiat 

cosmology and iJo=100 km s~^ Mpc^^, (7o=0.5. (j){M,z) 
is the luminosity function in the Schechter (1976) 
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Fig. 4. — Change in (UV— B) color vs. redshift, normalized to 
2 = 0. The calculations are based on the HD 165763 WC and the 
WR152 WN stellar spectra. 
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parameterization, 0(M) = OAlnlOcj)* [iqO a{m' -M)^{i+a) 
ea;p(— 10°'^^^ ~^^), where (f>* is the normalization (galax- 
ies per unit volnmc), a. is the faint end slope, and M* 
is the characteristic absohite magnitude. To convert the 
SSRS2 Sm-Irr M* = -19.78 to the UV bandpass, we add 
the median (UV— i?)o = —0.91 color for the UV galaxies. 
To simulate WR galaxies, we add the relative change in 
(UY-B) over redshift (Figure 4): 

M;*jv = M*s + (UV - B)o + A(UV - B) [z] (2) 
To find the number of galaxies in redshift bin {z,z + 6z), 
we integrate 4'{M, z) to the limiting absolute magnitude 
Miim{z) over the volume of the redshift interval z to z+6z. 
The selection functions represent the number of galaxies 
expected per unit redshift if galaxies are uniformly dis- 
tributed. 

The curves in Figure 2 are the UV selection functions 
normalized to the observed number of galaxies. The selec- 
tion function for the S2000 LF (M* = -20.59, a = -1.51) 
is very similar to the SSRS2 Sm-Irr LF {M* = -19.78 + 
(UV-5)o, a = -1.81), and matches the (UV-5)o > -2 
galaxy redshift distribution (Figure 2, upper panel) quite 
well. 

In contrast, the UV selection function we calculate for 
WR galaxies predicts more WR galaxies at large redshifts. 
The redshift distribution of (UV— i?)o < — 2 galaxies (Fig- 
ure 2, lower panel) matches the WR selection function 
curves much better than the S2000 curve. The shifts we 
make in M* reflect only the change in the flux through 
the UV bandpass as a function of redshift. Thus the UV 
emission lines of WR galaxies can explain the extreme 
blue (UV— B)o colors and the redshift distribution of the 
(UV-B)o < -2 galaxies. 

We calculate the number density of Century Survey 
(which includes the SA57 region) galaxies to the same ef- 
fective limiting magnitude and solid angle as the UV sam- 
ple, and find that the extreme blue objects are ~1% of a 
"normal" (i?-selected) galaxy population. It is difficult to 
compare the known WR galaxy fraction. We estimate a 
lower limit by assuming the WR catalog of Schaerer et al. 
(1999b) is complete to i?=14 mag over the northern sky 



above 15 deg galactic latitude. The resulting WR galaxy 
fraction lower limit is 0.75%. Thus the extreme blue galaxy 
number density is consistent with the abundance of known 
WR galaxies. 

5. CONCLUSIONS 

We show that the extreme blue (UV— B)o < —2 galaxies 
observed by S2000 may be WR galaxies, starburst galaxies 
with strong UV emission from WR stars. The model ex- 
plains the observed (UV— S)o < —2 colors and the redshift 
distribution for the bluest objects. 

The average {V — R} color of 13 known WR galaxies 
taken from the literature, (V — R) = 0.27 ± 0.34, is con- 
sistent with the range of (V — R) colors we observe for the 
extreme blue (UV— B)o < —2 galaxies. We measure a me- 
dian optical color (V-R) = 0.38 ± 0.06 for 67% of the 
S2000 SA57, UV-selected sample. 

The UV emission lines can double or triple the flux 
through the S2000 UV bandpass at large redshift. As a re- 
sult, the (UV— B) color of a WR galaxy can become up to 
1.3 mag bluer at larger redshift. This changing (UV— -B) as 
a function of redshift predicts a selection function weighted 
toward large redshifts. Not only do we find an excess of 
(UV— i?)o < —2 galaxies at large redshift, compared to 
the (UV— i?)o > —2 sample, but our model matches the 
observed redshift distribution well. We plan to test this 
model by obtaining high quality spectra of the suspected 
WR galaxies and looking for the telltale broad WR emis- 
sion features. 
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